First observation of CP violation and improved measurement of the branching 
fraction and polarization of B° — > D* + D*~ decays 
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We report the measurement of the branching fraction, the polarization, and the parameters of the 
time-dependent CP violation in B° — > D* + D*~ decays using a data sample of 772 x 10 6 BB pairs, 
collected at the T(4S') resonance with the Belle detector at the KEKB asymmetric-energy e + e~ 
collider. We obtain a branching fraction of B = (7.82 ± 0.38 ± 0.63) x 10~ 4 , a CP-odd fraction 
of R± — 0.138 ± 0.024 ± 0.006 and, additionally, a fraction of the longitudinal component in the 
transversity base of _Ro = 0.624 ± 0.029 ± 0.011. The measured values of the parameters of the CP 
violation are S D *+ D *- = -0.79 ± 0.13 ± 0.03 and A D *+ D *- = 0.15 ± 0.08 ± 0.04. 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 



The analysis of the time-dependent decay rate of neu- 
tral B mesons allows for the measurement of CP vio- 
lation in the B system, which is related to the complex 
phase in the Cabibbo-Kobayashi-Maskawa (CKM) quark 
mixing: matrix The decays of a neutral B meson 

to two charged charm mesons have been analyzed in nu- 
merous final states both by the Belle and the BaBar col- 
laborations The mixing- induced CP violation in 
these decays is related to one of the angles of the CKM 
triangle fa = arg {-V c dV* b /V t dV£), which was precisely 
measured in b — > (cc)s decays [H y|- 

The time-dependent decay rate of a neutral B me- 
son, originating from the decay of an T(4S'), to a CP- 
eigenstate is given by 

P B o (At) = - e /t b'> 

4t b o 

x {l + <7[S'sin(AmAi) + ^cos(AmAt)]}, (1) 

where q is either +1 or —1 when the non-signal B meson 
decays as B° or B° respectively, At is the decay time 
difference of the two B mesons, tqo is the lifetime of the 
neutral B meson, and Am is the mass difference of the 
two mass-eigenstates of the neutral B system. 

The decay B° -> D* + D*~ is the decay of a pseudo- 
scalar to two vector mesons and the final state is not a 
pure CP-eigenstate but a mixture of CP-even and CP- 
odd final states, depending on the relative angular mo- 
mentum of the D* mesons Q. For an angular momen- 
tum of zero or two, the final state is CP-even while, for 
an angular momentum of one, it is CP-odd. The dom- 
inant contribution to the decay B° —> D* + D*~ is the 



tree-level b — > ccd transition. Contributions from pen- 
guin diagrams are possible but strongly suppressed in 
the Standard Model 0. 

The measurement is performed using the final Belle 
data sample consisting of an integrated luminosity of 
711 ftT 1 containing (772 ± 11) x 10 6 BB pairs, collected 
at the T(4S) resonance at the KEKB asymmetric-energy 
e + e~ collider The electron beam has an energy of 
8 GeV and the positron beam 3.5 GeV, which leads to a 
boost of the T(4S*) of /?7 = 0.425 along the beam axis. 
Therefore, the neutral B mesons have an average abso- 
lute flight length difference of about 200 fim, which allows 
for the determination of the decay time difference At. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of aero- 
gel threshold Cherenkov counters (ACC), a barrel- like 
arrangement of time-of-flight scintillation counters, and 
an electromagnetic calorimeter composed of CsI(Tl) crys- 
tals (ECL) located inside a superconducting solenoid coil 
that provides a 1.5 T magnetic field. An iron flux- return 
located outside the coil is instrumented to detect K\\ 
mesons and to identify muons. The detector is described 
in detail elsewhere jllj. Two inner detector configu- 
rations were used. A 2.0 cm beam-pipe and a 3-layer 
silicon vertex detector were used for the first sample of 
152 x 10 6 BB pairs, while a 1.5 cm beam-pipe, a 4-layer 
silicon detector and a small-cell inner drift chamber were 
used to record the remaining 620 x 10 6 BB pairs [l2|]. In 
2011, the data set recorded with the second configura- 
tion of the SVD was reprocessed using new track finding 
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algorithms. This significantly improved the track recon- 
struction efficiency, especially for low momentum tracks. 
Due to the presence of up to two charged slow pions in 
the final state and the high multiplicity of the decay, this 
leads to an increase in the reconstruction efficiency of 
B° — > D* + D*~ decays of 79% in comparison to the last 
Belle measurement, which used 657 x 10 6 BB pairs [5j. 

The analysis of B° — > D* + D*~ decays is performed 
by reconstructing charged D* mesons via D* + — > D°tt + 
and D* + — > D + tt°. Charged D mesons are recon- 
structed in the modes K~tt + tt + , K s tt + , K s ir + ir°, and 
K + K~ir + , and neutral D mesons in K~ir + , K~ir + ir°, 
~ and K + K~. Neutral pions are 
reconstructed in the ir° — > 77 decay mode from photons 
with an energy greater than 30 MeV and are required 
to have an invariant mass that lies within 15 MeV/ c 2 of 
the nominal mass [13J; this corresponds to a width of 



3.3 a . The neutral kaons, reconstructed in K s — > ir + ir 
decays, must have an invariant mass in the range 470 
MeV/c 2 < < 520 MeV/c 2 and fulfill the criteria 

described in Ref. [lH • Charged tracks used for the recon- 
struction of D mesons are required to have a distance to 
the interaction point along (perpendicular to) the beam 
direction of less than 4 (2) cm. We separate charged 
kaons and pions based on the information from the ACC, 
the time-of- flight measurement, and the measurement of 
the energy loss in the CDC. The applied selection has 
an efficiency of 98% (97%) and a misidentification rate 
of 9% (18%) for pions (kaons). Additionally, an electron 
veto based on the measurement of the shower shape and 
the energy deposit in the ECL is applied. Pions that are 
used to form D* mesons are denoted as slow pions, due 
to their low momentum, and need not pass particle iden- 
tification and the distance requirement because of their 
low spatial resolution and short range in the detector. 

Final states containing more than one D* + — > D + tt° 
candidate or more than one Kg candidate are discarded 
due to their high background level. D and D* candidates 
are selected based on their invariant mass and the mass 
difference between the D* and the D meson candidate. 
The applied selection is decay mode dependent and has 
an efficiency of about 95% for the decays of D mesons to 
two charged tracks and 85 to 90% otherwise. These re- 
quirements as well as the selection criteria for kaons and 
pions are determined using Monte Carlo (MC) simulated 
events by optimizing a figure-of- merit S/\/S + B, where 
S is the number of signal events and B is the number 
of background events. Charged tracks originating from 
a D meson decay are constrained to have a common ori- 
gin. The slow pions are constrained to originate from 
the decay vertex of the D* mesons, which is obtained by 
projecting the D momentum vectors to the beam inter- 
section region. Two D* candidates are combined to form 
a B° candidate. Its decay vertex is determined by a kine- 
matic fit of the tracks of the two D mesons to a common 
vertex with the constraint that they originate from the 
beam interaction region. 

Two variables, the beam-constrained-mass M^ c = 



1/4 — (p*) and the energy difference AE = E* B — 

with being the beam energy in the center-of-mass 
system, p£ the momentum, and E* B the energy of the 
fully reconstructed B meson, also in the center-of-mass 
system, are used to discriminate between signal decays 
and the background. Events with M^ c < 5.23 GeV/c 2 or 
|A_E| > 0.14 GeV are rejected. We also define a tighter 
selection region, referred to as the signal region, with 
A/bc > 5.27 GeV/c 2 and \AE\ < 0.04 GeV. This region is 
used to optimize the selection criteria as described above. 
In 32% of the signal events, there are multiple candidates. 
In the case of multiple candidates, the best candidate is 
selected based on the masses of the D mesons and the 
mass differences between the D* and D meson candi- 
dates. 

The selected sample consists of combinatorial back- 
ground and signal. Not all signal events are com- 
pletely correctly reconstructed: in some, a neutral pion is 
wrongly reconstructed or a charged kaon is misidentified. 
This component, denoted as cross-feed, is treated as sig- 
nal and fitted with its own probability density function 
(PDF). 

The number of reconstructed signal events, used for the 
calculation of the branching fraction, is determined by 
an extended, two-dimensional, unbinned maximum like- 
lihood fit of the Mbc and AE distributions. The signal is 
described in AE by the sum of two Gaussians and a bi- 
furcated Gaussian with a common mean, and in M\, c by 
an empirically determined parameterized signal shape in- 
troduced by the Crystal Ball collaboration [15]. The frac- 
tions of the individual terms and the ratio of the widths 
are determined using simulated signal events, while the 
mean and a scale factor of the width, which is introduced 
to absorb the possible difference between data and sim- 
ulation, are floated in the fit. The cross-feed is described 
by an additional term consisting of the sum of a con- 
stant term and a Crystal Ball function in AE, and a 
single Crystal Ball function in M^ c . Its shape and frac- 
tion are determined using simulated signal events and 
fixed in the fit. The fraction of the cross-feed is found 
to be 11.6% of the total signal yield. The combinato- 
rial background is described by a second-order polyno- 
mial in AE and an empirically determined parametrized 
background shape introduced by the ARGUS collabora- 
tion [16[ in Mbc- The fitted number of signal events is 
1225 ± 59. The fit projections and data distributions 
are shown in Fig. [TJ Together with the reconstruction 
efficiency, obtained using Monte Carlo simulated signal 
events, and the branching fractions of the decays of D* 
and D mesons, taken from Ref. [13j, this gives a branch- 
ing fraction of B(B° ->■ D*+D*~) = (7.82 ± 0.38) x 10" 4 . 

The contributions to the systematic error of the 
branching fraction are listed in Table [I] The dominant 
errors are the uncertainty of the reconstruction efficiency 
of charged tracks, neutral pions, and slow pions, and 
of the selection performed for the separation of charged 
kaons and pions (PID). These uncertainties have been 
estimated using high-statistics control samples. The fit 
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TABLE I. Systematic errors of the branching fraction. 
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FIG. 1. Fit projections and data points of Mbc (top) and AE 
(bottom). The solid line shows the total fitted distribution in 
the signal region of the other variable. The dashed line shows 
the fitted background distribution. 



model error is determined by varying all parameters that 
are obtained from simulation and fixed in the fit by two 
standard deviations (<r) to account for the statistical un- 
certainties of the Monte Carlo samples and possible sys- 
tematic differences between data and simulation. We re- 
peat the fit for each parameter and add all deviations 
from the nominal fit result in quadrature, where the 
fraction of the cross-feed is varied by 20% of its value. 
The event reconstruction efficiency error includes the un- 
certainty of reconstruction efficiencies as well as uncer- 
tainties in efficiencies of selection criteria related to the 
masses of the D and D* mesons. The uncertainty in 
polarization has a negligible effect upon the overall sys- 
tematic error. The total error is calculated by adding all 
contributions in quadrature. 

In total, we measure a branching fraction of 
B( B °^d*+d*-) = (7.82 ±0.38 ±0.63) x 10~ 4 , where the 
first error is statistical and the second is systematic. This 
is consistent with previous measurements made by the 
Belle, BaBar, and CLEO collaborations, which were per- 
formed using data samples having sizes of 1 — 30% the 
size of our sample [17H19l | . 

For the measurement of the parameters of CP viola- 



Source 


Systematic error (%) 


Charged track reconstruction 


±1.7 


jKs reconstruction 


±0.8 


7T° reconstruction 


±3.0 


Slow pion reconstruction 


±3.2 


PID selection efficiency 


±5.0 


N B b 


±1.4 


Fit model 


±2.1 


D and D* branching fractions 


±3.1 


Event reconstruction 


±0.8 


Total 


±8.1 



tion, additional requirements are made as described in 
Ref. @. The parameters S and A are determined with 
a simultaneous, five-dimensional fit. In addition to M^ c , 
AE, and At, we fit two angles in the transversity basis 
to measure the polarization and to statistically separate 
the CP eigenstates. We define the transversity basis with 
the rc-axis pointing in the direction of the momentum of 
the D* + meson in the CMS, and the decay products of 
the D*~ meson lying in the xy-plane in the rest frame 
of the D* + meson. The angle tr is defined as the angle 
between the momentum of the pion from the D* + decay 
and the z-axis in the D* + rest frame, and 0\ as the angle 
between the momentum of the pion from the D*~ decay 
and the x-axis in the rest frame of the D*~. The angular 
distribution of the B° decay products is given by 



1 d 2 T(B° -> D*+D* 

T G?COS0 tr <icOS0i 



i=0,X,i 



where 



H = 2 cos 6*i sin tr , 
Hx = 2 sin 2 6»i cos 2 6> tr , 
H u = sin 2 6>i sin 2 tr 



COS0tr, COSfij, 

(2) 

(3) 
(4) 
(5) 



and Rq, Rj_, and R„ are real- valued parameters that are 
defined to have a sum of one [2(|. The term containing 
R± represents the CP-odd component, while the other 
two represent the CP-even components. The measured 
distributions differ from this theoretical expectation due 
to mis-reconstruction, the angular resolution of the slow 
pions, and varying reconstruction efficiency. The shapes 
of the three components are obtained separately from 
simulated signal events. The signal distributions are de- 
scribed by fourth-order polynomials, where the terms of 
odd order are fixed to zero in cos0t r . The parameters Ri 
are corrected for the relative reconstruction efficiencies 
of the different components. The background in cos0t r is 
described by a polynomial of second order with the first 
order fixed to zero, and in cos 9\ with a polynomial of 
fourth order. 

The decay vertex of the accompanying B meson (re- 
ferred to as the tag B) is obtained from charged tracks 
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that are not used for the reconstruction of the signal B 
meson. This is done by constraining them with a kine- 
matic fit to originate from a common vertex and discard- 
ing tracks that likely originate from secondary decays. 
The flight length difference is given by Az = zcp — zt a ,». 
The procedure for flavor tagging is described in Ref. [21j . 
For each candidate, it gives the flavor q of the tag B me- 
son and a tagging quality variable r ranging from r = 
for no flavor discrimination to r = 1 for unambiguous fla- 
vor assignment. The event sample is divided into seven 
bins of r. The wrong tag fraction w and the wrong tag 
fraction difference between the tagging of B° and B° 
mesons Aw were determined in each of the bins using 
high statistics control samples @ . Equation Q] is modi- 
fied to account for the influence of imperfect flavor tag- 
ging and the dilution by different CP-modes in the final 
state to 



V B o (At,cos8 tI ,cos6i] 



1 



Atro 



I 



x (1 - qAw + q(l - 2w){[l - 2P odd (cos0 tr , cos0i)] 

x Ssm(AmAt) + A cos(AmAt)}) . (6) 

This modification includes the assumption that S+ = 
—S- = S and A + = A- = A, where S+ and A + describe 
the CP violation in the CP-even and S- and A_ in the 
CP-odd component. P dd is the probability of an event 
to be CP-odd, derived from the angular distributions and 
the parameters Ri with 

P_LP±(c0s6>tr,C0s6»i) 
Podd(C0S6 , tr,C0S6'i) = — — 

2 , R l H i (cos9tr,cos9i / 

i=0,_L,ll 



(7) 



The values of Am and r^o are fixed to the current world 
averages of 507 ps -1 and 1.525 ps, respectively This 
signal PDF is convolved with a term describing the de- 
tector resolution and effects from secondary decays of the 
tag B. The resolution function and its parameters were 
obtained from high statistics control samples and are de- 
scribed in detail in Ref. [6]. We multiply the resulting 
PDF with Eq. [2] and the according PDF in M hc and AE 
to obtain the overal signal PDF. The background in At 
is modeled with the sum of an exponential and a prompt 
component, convolved with the sum of two Gaussian dis- 
tributions whose widths depend on the uncertainty of 
Az. 

The fraction of signal events is determined separately 
for each bin in r using the fit model in AE and M\, c , 
which is the same as that used in the branching frac- 
tion measurement. In total, the PDF describing the 
signal and background in My,c, AE, cos6* tr , cos6*i, and 
At has 30 free parameters, with four of them being 
S, A, Rq, and R±. They are determined with an un- 
binned maximum likelihood fit. We find Po = 0.624 ± 
0.029, R± = 0.138 ± 0.024, S = -0.79 ± 0.13, and A = 
0.15 ± 0.08, with the statistical correlations given in Ta- 
ble [TTJ Figure [5] shows the fit projections and the data 



points for the two angles. Figure [3] shows the At dis- 
tribution of well-tagged events (r > 0.5) with q = 1 
and q = —1 and the corresponding raw asymmetry 
(N+(At) - N-(At)) I (N+(At) + AT_(Ai)), with Ni(At) 
being the number of events with flavor tag i in decay time 
bin At. 



TABLE II. Statistical correlations of the physical parameters. 
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FIG. 2. Fit projections and data points of cos# tr (top) and 
cos 0i (bottom) in the signal region of Mb c and AE. The 
black, solid line shows the total fitted function, the green, 
long dashed line the fitted signal, and the red, short dashed 
line the background. The brown dotted, the blue dash dotted, 
and the cyan dash double dotted lines show the contributions 
of the Ro, R±, and R, components, respectively. 

The systematic errors are summarized in Table llTfl Ma- 
jor contributions come from the signal resolution func- 
tion, the vertex reconstruction, and tag side interference. 
The errors due to the uncertainties of the vertex res- 
olution function are studied by varying each fixed pa- 
rameter obtained from data (simulated events) by one 
(two) standard deviations and repeating the fit. The 
interference of Cabibbo-favored b — > cud and doubly 
Cabibbo-suppressed b — > tied decays of the tag B meson 
is noted as tag side interference and described in detail 
in Ref. [22| . Its influence is estimated by generating ran- 
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TABLE III. Systematic errors of S, A, R , and R± 



FIG. 3. The top plot shows events with r > 0.5 for q = 1 and 
g = — 1 in the signal region of and A_E. The red, dashed 
(blue, solid) line and markers show the data points and the 
corresponding fit of events for q = 1 (q = —1). The thin, 
black line shows the estimated background contribution. The 
bottom plot shows the raw asymmetry of the two top curves. 



dom data samples from a modified PDF that includes 
tag side interference and fitting them with the nominal 
PDF. The difference between the fit result and the in- 
put value is taken as the systematic error. The vertexing 
error combines multiple sources of error such as the re- 
quirement for the vertex fitting quality variable, the se- 
lection in At, and the selection criteria for the tracks used 
for the determination of the decay vertex of the tag B. 
It also contains systematic errors due to imperfect SVD 
alignment and potential bias in the measurement of Az, 
which are estimated with Monte Carlo simulations. The 
fit model error has been estimated in the same manner 
as done for the branching fraction measurement. The 
contributions due to the physics parameters, the flavor 
tagging model, and the reconstruction efficiencies of the 
three polarizations are all estimated by varying the cor- 
responding values within their uncertainties. The total 
systematic error is calculated by adding all contributions 
in quadrature. The total systematic error on S is re- 
duced significantly in comparison to the previous Belle 
measurement by including R± as a free parameter in the 
fit for S and A. 

We estimate the significance of the fit according to 
Wilks' theorem [23]. The hypothesis of no CP violation 
(S = A = 0) is excluded with a significance of 5.4 a. This 
significance takes into account the systematic uncertainty 
by convolving the likelihood function with a Gaussian 



S 



A 



Ro 



Ri 



Fit model 

Physics parameters 

Flavor tagging 

Tag side interference 

At signal resolution 



±0.002 < 0.001 
±0.004 ±0.001 
±0.003 ±0.002 
±0.007 ±0.032 
±0.021 ±0.006 



Reconstruction efficiencies < 0.001 < 0.001 
Vertexing ±0.017 ±0.021 



±0.010 
±0.001 
< 0.001 
±0.002 
±0.001 
±0.002 
±0.004 



±0.003 

< 0.001 

< 0.001 
±0.001 
±0.001 
±0.001 
±0.004 



Total 



±0.029 ±0.038 ±0.011 ±0.006 



distribution whose width is equal to the total additive 
systematic error. The measurement is consistent with the 
Standard Model expectation of A w and S = sin 2</>i @ 
within 1.4 a. 

Additionally, we release the assumption, that S+ = 
— S- and A + = A-. This allows for the possibility of 
different relative penguin contributions to the two differ- 
ent components in the fit. The according signal PDF in 
At is given by 



r B o(At) 



1 -|At|/ 
p > t r0 



b° (1 - qAw + q(l - 2w) 



4t b o 

x {[(1 - P od d)S+ ± PoddS-} sin(AmAf) 

± [(1 - Podd)A+ + P ddA-} cos(AmAi)}). (8) 

Note the difference in the sign convention of SL in com- 
parison to Sj_ in Ref. 0. We obtain S+ = -0.81 ±0.13 ± 
0.03, A+ = 0.18 ± 0.10 ± 0.05, S_ = 1.52 ± 0.62 ± 0.12, 
and A_ = -0.05 ± 0.39 ± 0.08. Within the uncertainties 
we find no indication for different CP violation in the 
CP-even and CP-odd component. 

In summary, we report the measurement of the branch- 
ing fraction, the polarization, and the parameters of time- 
dependent CP violation of B° -> D* + D*~ decays. The 
results are 

B (B »^D*+D>-) = [7.82±0.38(stat.)±0.63(syst.)] x 10~ 4 
P = 0.624 ± 0.029(stat.) ± 0.011(syst.), 
R x = 0.138 ± 0.024(stat.) ± 0.006(syst.), 
S = -0.79 ± 0.13(stat.) ± 0.03(syst.), and 
A = 0.15 ± 0.08(stat.) ± 0.04(syst.). 

These are consistent with previous measurements by the 
Belle collaboration 0, [llj and supersede them. It is also 
in agreement with the Standard Model expectation and 
is the first measurement of B° — s- D* + D*~ decays that 
exhibits CP violation with a significance greater than 
5cr. 
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